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Abstract 

Building on the results of Ma, Trudinger and Wang [T7] , and of the author [16] , we 
study two problems of optimal transportation on the sphere: the first corresponds to the 
cost function cP{x,y), where d{-, •) is the Riemannian distance of the round sphere; the 
second corresponds to the cost function — log \x — y\, it is known as the reflector antenna 
problem. We show that in both cases, the cost-sectional curvature is uniformly positive, 
and establish the geometrical properties so that the results of [T6] and [T7] can apply: 
global smooth solutions exist for arbitrary smooth positive data and optimal maps are 
Holder continuous under weak assumptions on the data. 
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1 Introduction 



1.1 Monge-Kantorovitch problem on a Riemannian manifold 

Let M be a topological space, let c : M x M — )■ M U {+00} be a cost function and yUo^/^i 
be probability measures on M. In the optimal transportation problem, one looks for a map 
T : M — )■ M that minimizes the functional 

(1) ^(T)= / c(x,T(x))rf/io(x), 

J M 

under the constraint that T pushes forward /iq onto /^i, (hereafter T^/io = /^i), i.e. 

MB dM Borel, = /io(T"^(5)). 

This problem has been first studied by Monge [in], with M = M", for the cost c = |x — ?/| 
(the Euclidean distance). Beyond Monge's original problem, Brenier studied the case of the 
quadratic cost c = |x — yp, and pointed out its close connection with important nonlinear PDEs 
(Monge-Ampere, Euler etc.). For the quadratic cost, when /io is absolutely continuous with 
respect to the Lebesgue measure, he proved the existence and uniqueness of an optimal map 
T. This maps has a convex potential (i.e. T = V0 with convex) and is shown to be the only 
map with convex potential that pushes forward /zq onto /ii. After Brenier 's result, the theory of 
optimal transportation has been extended to general cost functions. The existence of optimal 
maps is granted under very generic conditions on the cost function, and the way to obtain 
it is achieved through a general procedure, known as Kantorovitch duality: Optimal maps 
are obtained by solving the dual Monge-Kantorovitch problem, whose unknown are potential 
functions. For (p a lower semi-continuous function on M, we define its c-transform as 

4>^{x) = sup{— c(x, y) — 0(x)}. 

j/gM 

A potential (p is c-convex if it is the c-transform of some other ip : M ^ M.. In that case, the 
equality (j) = (p^^ holds. (Notice that the quadratic cost is equivalent to the cost —x ■ y, for 
which the c-transform is nothing but the Legendre-Fenchel transform; hence, [—x ■ yj-convex 
functions are convex functions.) Under suitable assumptions, and following for example [B], the 
minimizers in ([1]) are related to c-convex potentials as follows: for an optimal Topt in ([1]), there 
exists a c-convex potential (p such that 

(2) for a.e. x G M, T^ptix) = G^{x) := {y € M, (P{x) + <j)\y) = -c{x, y)}. 

Conversely, if T : M — )■ M can be expressed under the form ([2]) for some c-convex (p, for /xq a 
probability measure on M and fii its push-forward by T, then T is the optimal map between 
^0 and /ii. Of course, it is not clear a-priori that ([2]) defines a map, as G(f,{x) is a set. However, 
under a suitable assumption on the cost (assumption Al below), the set G^{x) will be reduced 
to a single point for Lebesgue almost every x. Note also that when M is compact (which we 
will assume throughout the remainder of the paper), the set G^{x) is never empty when (p is 
c-convex. 
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Brenier's result was generalized in a natural way to Riemannian manifolds by McCann: Let 
M be a manifold, with Riemannian metric g, compact and without boundary, with distance 
function d{-, •). For u,v E T^^M), {u,v)g{x) (or in short {u,v)g) denotes the scalar product on 
Tx{M) with respect to the metric g, = {v,v)g. From the results of [18], in the case where 
c = (P/2, the optimal map can be expressed as a so-called gradient map, i.e. 

G^x) = exp^(Vg0(x)), 

where Vp denotes the gradient with respect to the Riemannian metric on M (and from now, 
we omit the subscript g), and is some c-convex potential. For a general cost, one needs first 
to introduce the c-exponential map c-exp^(-), defined as the inverse of y — )■ — Va;c(a;,?/) (again 
see assumption Al below). Optimal maps will then be given by 

Topt{x) = G^{x) := c-exp^.(V0(x)). 

(This definition holds in the a.e. sense and is consistent with the definition ([2]).) In all cases, 
for a smooth potential such that G<^#/io = /ii, the conservation of mass is expressed in local 
coordinates by the Monge- Ampere type equation 

(3) det{D'<p + Dlc{^,G^)) = det D^,^c(x, 

where po, pi are densities of po, Pi with respect to the Lebesgue measure. 

A complete description of the optimal transportation problem can be found in [21] , otherwise 
the introduction of [16] encloses the necessary material for the present paper. We also mention 
the second book by Villani on optimal transport [25] , which presents some of the results enclosed 
in this paper. 

In the present work, we will address the problem of regularity of solutions of (or equiv- 
alently of minimizers of ([T])) in the particular geometrical setting of the constant curvature 
sphere of R"", that we will denote We will consider two cases: the quadratic cost c{x,y) = 

^(P{x, y) with d the Riemannian distance, and the reflector antenna case c(a;, y) = — log \x — y\ 
(see below). We will show that in both cases, the regularity results obtained in [TB] and [IT] 
hold. As noticed at the end of the paper, our result can be easily generalized to the case 
c{x,y) = f{d{x,y)), assuming some conditions on / (see Theorem 14.11) . 

The problem of regularity of optimal maps is related to the problem of regularity of solutions 
of the associated elliptic Monge-Ampere type equation. This fully non-linear elliptic partial 
differential equation has received considerable attention over the last decades. Up to recently, 
the only case where regularity results had been obtained were for the Monge-Ampere equation 

det D^u = /, 
u convex, 

with works by CaffareUi [D El [3l H El [7] , Urbas [23], and Delanoe [9]. This form of the Monge- 
Ampere equation is associated to the quadratic cost. It was only recently with the works of Ma, 
Trudinger and Wang [TTl [22] and subsequent results by the author ^JJ that regularity results 
for generic costs were obtained. Our goal in the present paper is to continue this study in the 
Riemannian setting, in the particular case of the round sphere. 
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The reflector antenna problem Consider a closed hypersurface S of parametrized by 
a so-called "height" function h from §"^^ to M, i.e. S = {xh{x),x G S"~^}. Given fi, fi' two 
domains of S"""^, and two probability measures fio, fJ'i on Q,Q', the reflector antenna problem 
is then to find h, under the constraint that the antenna reflects the incoming intensity po into 
the outgoing of intensity pi. If a ray with direction x is reflected into a ray of direction T{x), 
the conservation of energy turns into a Monge- Ampere type equation for h (see [8]): 

det {Vjjh + {u- ri)eij) _ po{x) 
r]"--^deteij pi{T{x)) 

with 

|V/lP + /l2 

Vijh the second covariant derivative of h, and Cij the Riemannian metric of S""^. 

Existence, regularity and uniqueness of weak and strong solutions have been addressed by 
several authors, among them Wang [26], Guan and Wang [H], Glimm and Oliker [13], and 
Oliker [20]. The connection with optimal transport theory was established by Wang in [27], 
who showed that the reflector antenna problem is equivalent to an optimal transport problem 
on the sphere with cost equal to — log \ x — y\. 

Organization of the paper In the next section, we will expose our results; then we will 
give a reminder of previous results and give some definitions and notations (some of them can 
be needed to understand the results). The rest of the paper is dedicated to the proofs of the 
results. 

2 Results 

This work addresses the issue of regularity of optimal maps. We deal only with global solutions 
(the measures and transport maps are defined on the whole of When the data are 

positive and (resp. C°°) smooth, we show that the optimal potential is (resp. C°°). 
When the target measure is bounded by below, and the source measure satisfies B(n-l), see 
([7j) (resp. A(n-l,p) for some p > n — 1, see (jS])) the optimal potential is (resp. C^'° for 
a{n,p)). We also give an original and self-contained proof of the connectedness of the contact 
set Gij){x) (defined in ([2])) in those particular cases. Our results are identical for the quadratic 
cost c(x, y) = ^(P{x, y) and the reflector antenna c(x, y) = — log \x — y\. 

For the reflector antenna problem, classical smooth solutions of (|3]) had been obtained in 
[2ni 123 [E] • Global regularity for weak solutions of (E]) had also been obtained independently 
in [S], under the assumption that both the source and target measures have densities bounded 
away from and infinity. Here we relax their assumptions, allowing the source measure to 
vanish, and requiring an integrability condition that does not even imply absolute continuity 
with respect to the Lebesgue measure. Moreover the Holder exponent in the (7^'°" result is 
explicit. 
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Our results are not a simple corollary of the results of [16] once the positivity of the cost- 
sectional curvature is asserted. Indeed, the basic assumption AO is not satisfied on x S""^, 
due to singularities of the cost functions (cut-locus), hence the equation itself could become 
singular. We have to show first that the graph of optimal transport map lies in a sub domain 
of §"~^ X S"^^ which is uniformly far from the cut locus. This is done by improving a result 
estabhshed in [TT], and also adapting it to the antenna case (see the Propositions 15.21 and 16. II) . 

Then we can localize the problem and reduce it locally to an Euclidean problem, and thus 
use the C^'" estimates of [IB] to show partial regularity under assumption ([H]) or ([7]). 



For classical regularity, we employ the method of continuity: combining the results of [IDj 

with the crucial a-priori estimate established in [T7|, we obtain classical smooth solutions. 

We choose to present in details the proof of the results for the quadratic case, from which 
the antenna case follows easily, once the key ingredients are verified. 

Acknowledgments I wish to thank Alessio Figalli and Cedric Villani for precious remarks 
and fruitful discussions. I also thank Philippe Delanoe with whom we started to think about the 
problem of regularity for optimal transportation on the sphere, and Robert McCann who first 
raised to me the issue of the connectedness of the contact set, in 2003. I gratefully acknowledge 
the support of a French Australian exchange grant PHC FAST EGIDE No.l2739WA. I wish also 
to thank Neil Trudinger, Xu-Jia Wang and the Center for Mathematics and its Applications at 
University of Canberra for their hospitality, as part of this work was accomplished there. 

2.1 The quadratic cost 

In this part, we consider S"~^ the unit sphere of M" equipped with the round metric g, and 
Riemannian distance d, and we let c{x,y) = ^cP{x,y). We first have the following essential 



Proposition 2.1 For all x G S"^"*^, the set {p G dxc{x,y),y G S"^"*^} is the closed ball 5(0, vr). 
Letting X be the antipodal point of x, the set {— Vxc(x, y), y G S"~^\x} is the open ball 3(0,7:) . 
Hence for all yo, yi G §"^^ \ x, the c-segment [?/o, yi\x (see definition \3.5\) is well defined. 

The next result shows that As is satisfied outside of the cut-locus. 

Theorem 2.2 Let antidiag be the set {{x,x),x G S"""*^} with x the antipodal point of x. Then 
the cost- sectional curvature is uniformly positive on S""-*^ x S"^-'^ \ antidiag. 

Then we show that the contact set (^^(x) is always connected (and even c-convex). The 
definitions of d(f){x),d^(j){x) are reminded below in Definitions 13.21 [3l3l 

Theorem 2.3 Let (p be a c-convex potential on S"^^ For all x G S"^"*^, let the set G^{x) be 
defined as in Then 

1. G^{x) is c-convex with respect to x, and d(f){x) = d'^(j){x). 

2. When x G G^{x), G^{x) is equal to the whole sphere §"^^. 




remark: 
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Finally, we conclude with our regularity result: 

Theorem 2.4 Let fio, fii be two probability measures on Ei"'^^ and cf) be a c-convex potential. 
Assume that G^^fio = yUi- Assume that fii > mdVol, for some m > 0. Then 

1. If fio satisfies B(n-l) (see then (p G C"'^(S"^"'^), and as in Theorem \3 . 1 1\ the modulus 
of continuity ofVcj) depends on f in 

2. If fiQ satisfies A(n-l,p) (see ^) for some p > n — 1, then (p G C^'^(S"^"'^) with [3 = 
P{n — l,p) as in Theorem VJ . 1 1\ 

3. ///io,/ii have positive C^'^ (resp. C°°) densities with respect to the Lebesgue measure, 
then (p G (:73."(§"-i) ^^^^y ^ ^ ^ ^^^^^^ ^ (:7°°(§"-i).; 



2.2 The reflector antenna 

The results presented above adapt with almost no modification to the refiector antenna. Before 
that, we remark that by changing y into — it is equivalent to study the cost c(x, y) = 
— log |x + ?/|. Then, as in the previous case, the set of singular points of the cost is equal to 
antidiag (i.e. the set of antipodal pair of points). Moreover, it is straightforward to check 
that whenever c{x,y) satisfies AO, Al, A2, AS, then c{x, —y) also satisfies those assumptions. 
Then, we prove the 



Theorem 2.5 Theorem \2. 3\ and \2.4\ hold for c{x, y) = — log \x — y 



(For point 2 of Theorem 12. 3^ notice that one must not consider x the antipodal point of x, but 
X itself.) Before entering into the proofs of our results, we present a review of previous results 
and concepts from the work of Ma, Trudinger and Wang in [Ej and the author in [1^. 



3 Cost-sectional curvature : from geometry of contact 
sets to elliptic regularity 

This paragraph is a short reminder of the results established in [TB], [22], [IZ]- In a nutshell, the 
things to remember are the following: when the contact set defined in ([2]) is unconditionally con- 
nected, the associated Monge- Ampere equation behaves well (i.e. smooth data lead to smooth 
solutions). Moreover, there is a tensor whose non-negativity is equivalent to the connectedness 
of the contact set. This tensor is the Ma- Trudinger- Wang tensor, or cost-sectional curvature 
tensor. When the cost is the squared Riemannian distance, this tensor contains informations 
about the Riemannian curvature of the underlying manifold. 

We adapt the results of [16], [22], [T7] to the Riemannian case by considering considering M 
a n-dimensional manifold, and D a general domain of M x M (which for simplicity we assume 
compact), instead of restricting to tensorial domains of x M", i.e. domains of the form ^2x^2' 
for fl, fl' domains of M". We denote tti, 712 the usual canonical projections. For any x G 7ii{D), 
we denote by the set D fl n^^{x). We proceed similarly for the y variable. Let us introduce 
the following conditions: 
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AO The cost function c belongs to 



Al For all x G 7Tx{D), the map y — ?■ — Vxc(x, y) is injective on D, 



A2 The cost function c satisfies det Dl 



,c 7^ for all (x, in 



Assumption Al allows to introduce the following definition: 

Definition 3.1 Under assumption Al, for x G 7ri(D) we define the a- exponential map at x, 
which we denote by c-exp^, the map such that 



The c-exponential map coincides with the Riemannian exponential map when c = d'^/2. In the 
general case, optimal maps are defined by 



for some c-convex potential 0. 

Noticing that for all y G G^{x), 0(-) + c{-,y) has a global minimum at x, we introduce / 
recall the following definitions: 

Definition 3.2 (subdifferential) For (p a semi-convex function, the subdifferential of (p atx, 
that we denote d(j){x), is the set 



The subdifferential is always a convex set, and is always non empty for a semi-convex function. 

Definition 3.3 (c-subdifferential) If (p is c-convex, the c-sub- differential of cp at x, that we 
denote d'^(p{x), is the set 



The inclusion ^ d^(p{x) C d(p{x) always holds. 

We make here the following important remark: c-exp^((9'^0(x)) is set of all y such that — 0(-) — 
c(-, y) reaches a global maximum at x, while c-exp^((90(x)) is the set of all y such that — 0(-) — 
c{-,y) has a critical point at x. 

We remind the definition of c-convexity (see [17]): 

Definition 3.4 (c-convex sets) Let x G 'Ki{D). A subset u' of tt2{Dx) is c-convex (resp. 
uniformly c-convex) with respect to x if the set {— Vxc(x, y),y E u'} is a convex (resp. uniformly 
convex) set of T^M. Whenever u x u' G D, u' is c-convex with respect to u if it is c-convex 
with respect to every x E uj. 



V(x,?/) G D, c-exp^{-Vxc{x,y)) = y. 



c-exp^(V0(x)) 



d(P{x) = {peT,{M) -.Wy eM,^{y)>(P{x) + {p-exp-\y)) + o{d{x,y))} 
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There exits also a notion of segment associated to the cost c, which we call c-segment: 

Definition 3.5 (c-segment) Let {x,yo) and {x,yi) G D. Let pi = —Vxc{x,yi),i = 0,1. 
Assume that [po,Pi] C —Vxc{x,Dx). Then the c-segment between yo and yi with respect to x, 
which we denote [yo,yi]x is defined by 

[yo,yi]x = {y e Dx,-Vxc{x,y) e [po,pi]}- 

Then yg will be such that —VxC^x^ye) = Opi + (1 — 0)pq. 

We now remind the definition of the cost- sectional curvature. Although this name and the 
present formulation comes from [16], this tensor has first been introduced by Ma, Trudinger 
and Wang in [17] . 

Definition 3.6 Under assumptions A0-A1-A2, for every {xo,yo) G D, one can define on 
TxgM X TxgM the real-valued map 



(4) &,{xo,yo){^, v) = {x,p) -> -c(x, c-exp^^^{p)) 



xo ,po=-'V xc{xo ,yo) 



When C,, y are unit orthogonal vectors (with respect to the metric g at xq), ©c(xo, yo){^, defines 
the cost- sectional curvature from xq to yo in directions (^, z^). This definition is intrinsic (i.e. 
coordinate independent) . 

We are now ready to introduce the following condition on the cost function: 

As (Positive sectional curvature) The cost-sectional curvature is uniformly positive on 
D. In other words, there exists Cq > such that for all {xq, yo) G D, for all ^,1^ E T^^^M, 



(5) e,{xo,yom,i^)>CoWm'. 

Alternatively, condition Aw (non-negative sectional curvature) is satisfied whenever 
(jSj) is satisfied with Cq = 0. 

In [H], [22], [ig the study was restricted to the case where M = R", D = x fi', for Q, Q' 
two domains of M". Concerning classical regularity. Ma, Trudinger and Wang [T7] proved the 
following: 

Theorem 3.7 (| I17] ) Let c be a C^{Q x Q') cost function that satisfies assumptions Al, A2, 
As on {Q X Q'), Q' being c-convex with respect to Vt. Let /io = podx,fii = pidx be probability 
measures respectively on Vt and Vt' . Assume that po,Pi are bounded away from 0, and belong 
to C'^{Q) (resp. C'^{Q')). Let (p be a c-convex potential on fl such that G^^^q = fii. Then 
(f) G for every a G [0, 1). 

The key step toward this result is an interior a-priori estimate for second derivatives of 
solutions of ([3]); we will rely on this a-priori estimate in the present paper. 

In [in] we proved the following relationship between regularity and a natural geometric 
property of the c-convex potentials: 
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Theorem 3.8 (| |16] ) Assume that for all fio,fii smooth positive probability measures, the c- 
convex potential (j) such that G^^fio = fii is . Then for all (p c-convex on Q, for all x ^ Q, 
d(j){x) = d'^4>{x). 

This resuh, combined with Theorem 13.71 yields easily the following: 

Theorem 3.9 (| |16] ) Let c be a C'^{VL x Vl') cost function that satisfies assumptions Al, A2, 
As on [Vt X Vt'), Vt^Vt' being c(resp c*)-convex with respect to each other. Then, for all 
c-convex on Q, for all x ^ Q, the set G^{x) is connected. 

In [16], it was also proved that the non-negativity of the cost-sectional curvature (i.e. con- 
dition Aw) is a necessary condition for regularity: without that condition, one can construct 
potentials that are not even though the data are smooth and positive. Relying on the results 
of [22j, it was then shown that the non-negativity of the cost-sectional curvature was equivalent 
to the connectedness of the set G^{x) (see ([2])). 

Theorem 3.10 (|]16], [22]) Let c be a C^(fi x Q') cost function that satisfies assumptions Al, 
A2 on {Q X Q'), Q,Q' being uniformly c(resp c*)-convex with respect to each other. Then the 
following are equivalent: 

1. The cost c satisfies Aw; 

2. For all c-convex, for all x &Vt, G(p{x) is connected; 

3. For all fiQ, /ii G°° -smooth positive probability measures on VL, VL' , the optimal map between 
yUo O'l^d Hi is 

Under positive cost-sectional curvature, we proved also C^'" regularity of optimal potentials 
for rough data. Let us first introduce some regularity properties on the measures /xqj/^i- The 
first one reads 

A(n, p) For some p Gjra, +oo], C^^ > 0, 

(6) ^^o{B,{x)) < C^oe"^^"^^ for all e > 0, a; G M. 
The second condition reads 

B(n) For some / : M+ ^ M+ with lim.^o /(e) = 0, 

(7) /io(5e(x)) < /(e)e"(^-") for all e > 0, x G M. 

We remark that h E L^ implies A(n, p) with the same p, while B(n) doest not imply fi E L^. 
Letting T-L^~^ be the n — 1 dimensional Hausdorff measure, we also notice that ([7]) implies 
that H does not give mass to sets A such that 'H"'~^{A) is finite, which is close to the optimal 
assumption for existence of an optimal map in ([T]) (e.g. /iq does not charge n — 1-rectifiable 
sets). We also denote classically Qs = {x E 0, : d{x, dQ) > 6). 
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Theorem 3.11 ([16]) Under the assumptions of Theorem \3.9{ on c, Q, Q' , let /io,yWi be prob- 
ability measures respectively on Q and Q' . Let cj) be a c- convex potential on Vt such that 
G^^f^fiQ = fii. Assume that fii > mdVol on Q' for some m > 0. 

1. Assume that fiQ satisfies A(n,p) for some p > n. Let a = 1 — ^, /3 = . Then for 
any 6 > we have 

U\\c^.p(ns) < C, 

and C depends only on 5 > 0, C^^ in on m, on the constants in conditions AO, Al, 
A2, As. 

2. If Ho satisfies B(n), then (j) belongs to C^iVts) and the modulus of continuity ofWcf) depends 
also on f in 

Moreover, it was shown that in the case where c is the squared Riemannian distance on M, 
denoting the sectional curvature of M, the identity 




holds (meaning that the two tensors coincide on (Tx^lY). Gathering all those results implies 
that for a manifold whose sectional curvature is negative at some point on some two-plan, one 
can exhibit positive smooth densities such that the optimal transport map between them (for 
the cost equal to d"^ {■,■)) is not continuous. 

4 Connectedness of the contact set: Proof of Theorem 

We prove the following more general result: 

Theorem 4.1 Let S"^^ be the unit sphere o/ equipped with the round metric g, and Rie- 
mannian distance d. Let c{x,y) = f{d{x,y)) for some / : [0, tt) — t- M smooth and strictly 
increasing, with /'(O) = and such that the conditions AO, Al, A2, As are satisfied on 
X gn-i ^ antzdzag. Then 

1. for all X, yo, yi G S""^ such that -x ^ {yo, yi}, for all y e [yo, yi]x, y 7^ -x; 

2. for all c-convex, for all y E §"^^, + c{-,y) reaches a global maximum at x = —y, and 
any other critical point is a global minimum; 

3. for all c-convex, for all x G S"^"*^, d(j){x) = cF(f){x). 
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Remark. We notice that d{x,y) = arccos(x ■ y), and that 

-\og\x + y\ = -i(log(^|x + yp) + cste 

= — - log(l + cos((i(x, y))) + cste. 

Hence in both cases the cost-function is of the form c{x,y) = f{d{x,y)) that fits into the 
framework of Theorem 14.11 

Proof of Theorem 14. 11 We first prove the point 1: as the cost c is under the form c = f{d), 
we have, for all x y, 



(8) V.c{x,y) = -f{d{x,y))e. 



yi 



where Cy = exp^^{y)/d{x,y). If the cost satisfies Al, then /' should obviously be a strictly 
monotone function. This yields that / is necessarily strictly convex. Now since / is non-negative 
and increasing, from the definition of the c-segment, we have easily using ([H]) 

max{d{x, ye), 9 G [0, 1]} < max{d{x, yo), d{x, yi)}- 

This yields point 1. Note that point 1 implies that the c-segment [yo, yi]x is well defined provided 
—X ^ {yo, yi}. We now prove point 2 and 3 for of the form 

(9) (j){x) = max{-c(x, yo) + Oq, -c(x, yi) + ai}. 

We let 

h{x) = + c(x, y), 

Dy = S--'\{-y}. 

Note that h is semi-convex on Dy. We say that x G Dy is a critical point of h whenever 
G dh{x). We now use the fact that the cost satisfies assumption As: 



Lemma 4.2 Under the assumptions of Theorem 4-1 • let x ^ —y he a critical point of h, then 
X is a local minimum of h. Moreover h reaches its global maximum at x = —y. 

Proof. We choose x ^ —y and assume that is not differentiable at x otherwise the conclusion 
is trivial. Without loss of generality, one can then write 

= max{-c(2;', yo) + c{x, yo), -c{x! , yi) + c{x, y^)}. 

Moreover, if + c{-,y) has a critical point at x where c{-,y) is differentiable, then necessarily 
y G [yo,yi]x, the c-segment with respect to x. Hence, one can write y = ye for some 9 G [0, 1]. 
Finally, if yo ^ yi, 4> is necessarily differentiable at x = —yi, i = 0, 1, since in a neighborhood of 
yo (resp. yi), = ~c{x',yi) +c{x,yi) (resp. = —c{x',yo) +c{x,yo)). Hence we can assume 
that X 7^ —yo and x ^ —y\. li y = yo or y = yi then (j) + c{-,y) has a global minimum at x. In 
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the case where y ^ {yo.yi} we use the result of [161 Proposition 5.1]. Locally around x, under 
assumption As, we have for G [e, 1 — e] 

max{-c(a;', yo) + c(x, ?/o), -c{x\ yi) + c{x, yi)} 
> -c{x', ye) + c{x, ye) + 5<f'{x, x') + o{(f{x, x')), 

where 5 > for e > 0, d{yi, y^) > 0. Hence, we have 

(j){x') + c{x', ye) > c{x, ye) + Scf{x, x') + o{(f{x, x')), 

with equality at x' = x. So if + c(-, y) has a critical point at x, either y = y^ or yi in which 
case X is a global minimum, oi y = ye for some 6 e (0, 1), and then x is a strict local minimum. 

This implies that h reaches its global maximum at x = —y: Indeed, if h reaches its maximum 
at X e Dy, then h is semi-convex around x (since it is the sum of a c-convex function and a 
smooth function), and x is thus a critical point of h. As we saw in Lemma \A.2\ any critical 
point of h in Dy is a local minimum, and as we are on a compact manifold, h must have a 
global maximum, which is thus necessarily —y. 

□ 

Lemma 4.3 Let h, (p be defined as above. Then any local minimum of h is a global minimum 
on E>''-\ 

Proof. We first assume that y ^ {yo,yi}- In this case, we prove the stronger assertion that 
h has only one local minimum (which in particular is the global minimum). So let us assume 
that h has two distinct local minima Xi ^ X2- Note that from Lemma [4.2[ both xi and X2 are 
different from —y, which is the global maximum of h. From the proof of Lemma [4. 2[ these two 
local minima are strict, and are thus separated: xi and X2 do not belong to the same connected 
component of a level set of h. Then consider 

inf sup h{-f{t)), 

where P is the set of continuous paths from [0,1] to §"^^ such that 7(0) = xi,7i = X2. 
By standard compactness arguments, this saddle point is attained at some point X3. Clearly 
^3 7^ ~y since, from Lemma S21 h reaches its global maximum at —y. Thus X3 G Dy, and X3 is 
a critical point of h. In view of Lemma W?2[ X3 should then be a local minimum of h, moreover, 
since y ^ {yo,yi}, this local minimum should be strict. We reach a contradiction. 
li y = yo, the situation is simpler: h reaches its global minimum on the whole set 

E = {x' : -c(x', yo) + > -c(x', yi) + aj. 

Moreover h reaches its global maximum at x' = —yi, and, from assumption Al, there are no 
critical points outside of E, so there can be no local minima of h outside of E. 

□ 

Lemma 4.4 Let be defined as above. Then for all x G S"~"^, d^(j){x) = d(j){x). 
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Proof. Consider p G d(t){x). Then p G [— Va;c(a;, yo), xc{x,yi)\x. If x is on the set where (f) 
is not differentiable, (otherwise there is nothing to prove) then —x ^ {yoiUi}- From the point 
1 already proved, there exists y G [yQ,yi\x C \ {— x}, such that p = —'Vxc{x,y). Then 
considering the function h = + c(-, y) on S"~^, clearly a; is a critical point of h. Thus, by what 
we proved above, a; is a global minimum of h, and so p G d'^cp^x). 

□ 

Then arguing as in [161 Proposition 2.11], if the conclusion of Lemma [4.41 holds for all c- 
convex functions of the form it holds for all c-convex; this achieves the proof of Theorem 

□ 

Remark. This argument allows to go from [TF, Proposition 5.1] to the important identity 
d^(j) = d(j). One of the crucial ingredients is the absence of boundary used in the saddle point 
argument. The case with boundary has been treated by Trudinger and Wang in [21]. An 
alternative argument to obtain d'^cf) = dcf) from Aw has been obtained in [15j . 



5 Proof of Theorem 2.4 



Strategy of the proof Most of the proof is contained in the following points: 

1- Given /io,/ii satisfying the assumptions of Theorem 12. 4[ there exists a constant a such 
that d{x,G^{x)) < tt — a for all x G S"~^. Hence, G(i){x) stays uniformly far away from the 
cut-locus of X. Then we can reduce locally the problem to an Euclidean problem. 

2- The assumption As is satisfied by the cost function distance squared on the sphere 

(Theorem [OD- 

Once this is established, we proceed as follows: 

Given xo G S"~^ we can build around xq a system of geodesic coordinates on the set 
{x,d{x,xo) < R} for R < n. From point 1-, for r small enough, the graph {{x , G ,f,{x)) , x G 
Br{xo)} is included in the set Br{xo) x -B,r-2r(a;o) on which the cost function is C°°. From point 
2- and using [T7], a C"^ smooth solution to on Br{xo) will enjoy a a priori estimate at 
Xo- This estimate will depend only on the smoothness of fj^o, fii, on r, and r is small but can 
be chosen once for all. 

Once a a priori estimate is established we use the result of [10] : the method of continuity 
allows to build smooth solutions for any smooth positive densities. 

Then, we use the results of [TB] to conclude our partial C^'" regularity result. 



5.1 Reduction of the problem to an Euclidean problem 

Uniform distance to the cut locus We show that there exists a subset S"^ of x S"^-*^ 
on which A0-A2 are satisfied, and such that the graph of G^, i.e. {{x,G^{x)),x G is 
contained in S^. This subset is defined by 

(10) Si = { (x, y) G §"-1 X §"-1 ■.dix,y)<n-a} 

where a > depends on some condition on fio,f^i- We first remark that 
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Lemma 5.1 For all x G Ei^^^ , the set {y, {x, y) G S'^} is c-convex with respect to x. 
Then we show the following crucial result: 

Proposition 5.2 Let fio,fii be two probability measures on let cf) be a c-convex potential 

such that G^^fiQ = /ii. Assume that there exists m > such that fii > mdVol and that fiQ 
satisfies B(n-l) (see Then there exists cr > depending on m and on f in such that 

{{x,Gif,{x)),x G §"^^} C Si, where is defined in /[IJ^) . 

Proof: We use pjj; in that paper, it was shown, for satisfying G^^/io = /ii, that we 
have \d(j)\ < vr — e, and e > depends on ||d/ii/dVol||Loo || {dfio/dVo\)~^ Considering 0"^ and 
G^o that pushes forward /ii onto /io, we also have a bound such as \d(j)\ < vr — e, where e > 
depends on ||d/io/dVol||2,oo|| (d/ii/dVol)~^ Here we slightly extend this bound to the case 
where /xq satisfies B(n-l) ([7]), and /ii > mdVol. 

The starting point of the proof is [111 Lemma 2] where it is shown that for all Xi,X2 G 
for all c-convex, 

d{G4x2),Xi) < 2n J , 

a[Xi, X2) 

where Xi is the antipodal point to Xi. Hence the set = {x : d{x,Xi) > 6} is sent by G^ in 
B^{xi) where 

d{G^{xi),Xi) 

€ = ZtT . 



Considering ip = (p^, we have G^#fii = /io- This implies 

Ho{B,{xi)) > HiiEs). 
Taking 6 > fixed (for example 6 = 7r/2), we have 

/^o(54<i(G^c(xi),£i)(a;i)) > Mni{D^), 
where is the half-sphere centered at z. Hence 

inf Un - d(f)%x)\] 

> inf |r : 3x G §""\ /io(54r(^)) > inf /ii(D^)} 
:= a. 

In particular, if /xq satisfies B(n-l), we have 

f^o{B,r{x)) < C(Vol(i?4.(a;)))^-^/("-^) 

for r small enough, hence if fii > mdVol, m > we find that a > and the uniform distance 
between G^{x) and the cut locus of x is asserted. Now notice that sup{\d(f){x)\,x G S"^"*^} = 
sup{\d(j)'^{x)\,x G S"~^}, and the proposition follows. 

□ 

Remark. It is enough to conclude the proof to assume that fio{B^{x)) — > uniformly with 
respect to e. 
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Construction of a local system of coordinates Given xq € S"~^, we consider a system of 
geodesic coordinates around Xq, i.e. given a system of orthonormal coordinates at xq and the 
induced system of coordinates on T^.q(§"'~^), we consider the mapping p G T^.q(§"'"^) — )■ exp^.^(p). 
This mapping is a diffeomorphism from -B/j(0) C M"~^ to Bb.{xo) as long as i? < tt. Then, for 
for r < cr/2 we have 

Si n {(x,?/) : X e Br{xo)} C Br{xo) X C 5^_2^. 

We now take r = cr/3 Hence, Br{0) is sent in i?,r-2r(0), and the cost function is C°° on 
Br{0) X _B^_2r(0)for r small enough. 

In this case, we have c-exp^(p) = exp^(p), and also 

(11) Dl^c{x,y) = -[D,exp.2-\y), 

where D^exp^ is the derivative with respect to f of f — )■ exp^(f). Assumption Al is trivial, 
since y is indeed uniquely defined hj y = exp^{p). From ffTTl) . assumption A2 is true on any 
smooth compact Riemannian manifold, since in this case Jac(f — )■ exp^ v) is bounded by above. 



5.2 Verification of assumption As: Proof of Theorem 12.21 

It has been established in [TT] that in a system of normal coordinates ei,..,en-i at x with 
ei = p/\p\, we have 

r iif2 = j = i 

I sm(|p|) 

This relies on the fact that Dl^c{x,exp^{p)) = [Dp{exp^{p))]^^Dx{exp^{p)), and follows by 
computations of Jacobi fields. Hence, we have, in this system of coordinates, 

T COS T 

M{x,p) := Dlc{x,c-exp^{p)) = /-(!- — ^)/', 

where r = \p\, I is the identity matrix of order n — 1, and I' is the identity matrix on 
vect(e2, .., e„_i). So, for a given u e M"~^, 

M(x,p) ■ v ■ V = l-yp - (1 - ^^"^^ )|Ht;P, 

smr 

where Hi; is the projection of w on p^. This can be written, using intrinsic notations 



(?^(c(,y))|,_p.(.).,.)„ = \v\l-{i- ' IT '^ )|n.i;, 





3 COS 


\z\ 




sin 
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where (■, ■)g denotes the Riemannian scalar product on c-exp^(§"~^) and = {p,p)g- 

Let a, /3 E R"'~\ t E M, ut = a + 1(3 and r = |ttf |. We assume that r < n, and Utv is the 
projection of w on m^. Assumption As is then equivalent to show 



^ [t ^ M(x, a + /3t) ■ ■ i;] < -Co|/3|>P, 
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for all a, ^ e c-exp^§" ^, v ± (3, which is equivalent to 



■^^ rcosr 1 ^ 

smr J 



for all a, ^ e ^, v J- /3. Without loss of generality, we assume by changing t in t — to that 
a J- /3. We have 

UtV = V - \Ut\~^{v ■ Ut)Ut, 

hence, using v ± (3, 

\Utv\^ = Ivl"" - 2{v ut)^\ut\-^ + (v ut)^\ut\-^ 

I |2 / \2| 1-2 

= \v\ — [v ■ a) \ut\ . 
Now, assuming 1^1 = 1 and {v, a)g = c, we have to evaluate ^F, where F = G o r, and 

G{r) =1 : )(1 - — , 

sm r 

r{t) = \a + tf3\. 

A computation shows that 

—F = — — < ■ n (r - smr cosr)(l -) 

r-^ I sm r 

1 2c^ 
H (sin r — r cos r) 

sm r r'* 

+ . o (smr - rcosr)(l -) 

I sm"^ r 

2 , 2c^ 
V — sin r cos r) 



sin^ r r^ 

1 . 6c2 
smr — rcosr 



smr ■ ■ r"^ 



Notice that we have for some C > 0, and for all r e [0, tt], 

(12) sinr — rcosr > Cr^, 

(13) r — sinr cos r > Cr^. 

Hence, all terms are positive, except the last one. The sum of the last two lines has for r e [0, tt] 
the sign of 

2(r^ — r sin r cos r) — 3(sin^ ^ — r sin r cos r) 
= ^ (r sin 2r + 3 cos 2r + 4r^ - 3) . 
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One can check (by two successive differentiations) that this quantity is non-negative for r G 
[0, ir]. Hence we have, using f[T^ [T5]) . 
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Remember that c = a ■ v, with l^;! = 1, hence |a| > |c| and r > |c|, to conclude that -^F > 
Co\f3\'. 

Hence in a normal system of coordinates at x, we see that As is satisfied at x for any y such 
that d{x,y) < vr. As seen in [16], the cost-sectional curvature is intrinsic. This proves Theorem 



□ 



5.3 C regularity for weak solutions 

Consider a system of geodesic coordinates around xq- As we have seen, using Proposition I5.2[ 
in this system of coordinates, for tq = cr/3, we have 

Moreover, c satisfies AO- As on -Bro(O) x i?^_2ro(0). We let = i?ro(0), R = n — 2ro, and 
Q' = Bji{0), and remark that Q,Q' are uniformly strictly c-convex with respect to each other 
if To is small enough. Now, given (p c-convex on S"~^ such that G^^fiQ = /^i, we consider (p the 
restriction of to Q. We recall the definition of the Monge- Ampere measure of (p with respect 
to /ii, defined by 

(14) Vfi Borel, G'f/ii(5) = /ii(G^(5)). 

We claim the following: 

Proposition 5.3 Let 0, 0, f2, fi' be as above. Let /xo,/^i be probability measures on such 
that G(f)^^Q = fii, and > mdVol; then 

1. (f) is c-convex on VL with respect to the restriction of c to VL x VL' ; 

2. coincides with G^ on Q; 

3. Gf dVol < idVol on n. 

(p — m 

Remark. Note that in general, G^^yUoln = /^iln' will not hold. Indeed, G^{Q) is a pri- 
ori strictly included il', hence for B C fi', /ii(-B) = Ho{G^^{B)) will be strictly larger than 
fio{G^^{B) n Q). However, point 2. suffices to conclude that = G^/ii. 

Proof. The first point is straightforward: indeed at every x & Q, (p has global c-support 
~4>^{y) ~ c(-,y), and y G il' since G^{Q,) C il'. 
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The delicate issue for the second part comes from the following point: we have 

G^ix) = I?/ : (p{x) + c{x, y) = mi{(f){x') + c{x' , y), x' G r^jj. 

As we see, the infimum above has to be on Q, whereas for the definition of G^p it has to be 
on §"^^. Thus the inclusion G(f){x) C G^{x) is straightforward while the reverse inclusion is 
unclear. To overcome this we use the Theorem 14.11 and the fact that for all (p c-convex on S^'^, 
d^(j) = d(f). Hence if + c{-,y) has a local minimum it has to be a global minimum in all of 
S"'"^. Hence, if y G G^{x), then + c{-,y) has a local minimum at x, and hence a global one, 
thus y G G(i,{x). 

Now the last point follows from the following lemma, which is a reformulation of [TSl Propo- 
sition 5.11] in the sphere case. 

Lemma 5.4 Let cj) be c-convex on S"^-*^. Assume that Gj/.^/io = /ii- Assume that /xi > mdVol 
on Q' . Then for all u C §"~^, we have 

f^oi^) > ^^ol{G^{u))), and hence, G^dVol < — Hq. 

Before proving this lemma, we conclude the proof of Proposition 15. 3^ by observing that, 
from the point 2 of Proposition 15. 3[ on fl we have G*dVol = G^dVol. 

□ 

Proof of Lemma [Ql In S""^ we consider N = {y e S""^ : 3xi X2 e E:'^~'^,G^{xi) = 
G^{x2) = y}. Then N = {y e S""^ : 0" is not differentiable at y}. Hence Vol(A^) = 0, 
and Vol^G^i^u) \ N) = Vo\{G^{u)). Moreover, on G<^(aj) \ A^, G'^^ is single valued. Then 
G^\G^{uj)\N) coo. Hence, 

/ioM > fioiG-^\G^iu)\N)) 
= fi,{G^{uj)\N) 
> mVo\{G^{u)\N) 
= mVo\{G^{uj)). 

We are now allowed to use [TBI Proposition 5.9] to prove Theorem 12.41 Let first recall the 
result: 

Proposition 5.5 (Proposition 5.9 [16]) Let Q,Q' be domains o/M"^^. Let be c-convex 
with G(i){Q) C Q' . Let c satisfy assumptions AO, Al, A2, As on Q,Q' with strictly 
c-convex with respect to each other, then 

• z/ G*dVol satisfies A(n-1, p) (see IjiS^), for some p > n — 1, then G Cl^^{Q), with 
I3{n — l,p) as in Theorem VJ . 1 1\ 

• z/GfdVol satisfies B(n-l) (seel^), then G 
This achieves the proof of regularity. 
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5.4 a priori estimates for smooth solutions 

Let be a C"^ smooth c-convex potential on S"^^ such that G^^fiQ = fii, with fio, having 
positive C^'^ smooth densities. We recall the a-priori estimate of [IZl Theorem 4.1], adapted to 
our notations: 

Theorem 5.6 ([17|) Let fi, fi' be subsets o/M". Let (p be a C"'(fi) c-convex solution of ^ 
with G^{Q) C Q' . Suppose assumptions AO, Al, A2, As are satisfied by c on Q x Q' . Then 
we have the a priori second order derivative estimate 

\D^(j)\{x) < C, 

where C depends on n, dist{x, dfl), sup^ 0, po, pi up to their second derivatives, the cost function 
c up to its fourth order derivatives, the constant Co in As, and a positive lower bound of 
I det Dl „c\ on fl. 

Applying this result with Q = Br{0), fi' = -B,r-2ro(0), we have a bound on D^(j){xo) that depends 
only on tq and on the bounds on po, pi. Finally, note that tq can be chosen once for all once a 
is known. 

Continuity method In [TU] it was established that given a C°° smooth c-convex potential 
0, a C°° smooth positive measure po, and pi = G^^/io, the operator 

F -.(f)-^ F{4>) = G^#po 

was locally invertible in C°° around pi. Then the existence, for a given pair of C°° smooth 
positive probability measures po, pi of a C°° smooth c-convex potential is granted once a- 
priori estimates for the second derivatives have been derived. Indeed, this follows from the 
concavity of the equation and and the well known continuity method (see [12]). 

We conclude the following: for po, pi having C^'^ smooth probability densities, there exists a 
(unique up to a constant) c-convex potential G C^'" for every a G [0, 1[ such that G^^/io = Pi- 
If moreover po, pi are C°°, G C°°. This concludes the third point of Theorem 12.41 

□ 

The proof of Theorem 12.41 is complete. 



6 Proof of Theorem 2.5 



For this case it has already been checked (see [26l \27\ [H]) that the cost function c{x,y) = 
— log \x — y\ satisfies As for x ^ y. We will just prove that under some assumptions on the 
measures po, pi we can guarantee that G^{x) stays away from x. This will imply that when 
the cost is — log |x -|- |/|, G^{x) stays away from x. Then, the proof of the Theorem 12.51 mimics 
the proof of Theorem 12. 4[ We prove the 
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Proposition 6.1 Let S" ^ be unit sphere o/M". Let c{x,y) = — log|x — y\. Let T : -"^ — t- 

§"^^ be a 2-monotone map, i.e. such that 

(15) \/xi,X2 e §"-\ c{xi,T{xi)) + c{x2, T{x2)) < c(x2, T(xi)) + c(a;i, T{x2)). 

Let /io, /ii be two probability measures on M". Assume that fj,i > mdVol for m > 0, and that /io 
satisfies B(n-l) Then there exists eo > depending only on m, f in ^ such that 

\/x e W'-\d{x,T{x)) > eo. 

Proof of Proposition 16.11 We follow the same lines as in [11]. From f lTSjl . we have 

log \xi - T{x2)\ < log |xi - T{xi) \ + log 2 - log \x2 - T{xi)\. 

Letting M = — log|xi — T(xi)|, the set {x : |x — T(xi)| > 1} is sent by T in the set {y : 
\y — xi\ < 2 exp(— M)}. Then the bounds on /iq, /ii yield an upper bound on M as in the proof 
of Proposition 15. 2[ 

□ 

Remark. The results that we have established for the reflector antenna and the quadratic 
cost can certainly be extended to the class of all costs that follow the assumptions of Theorem 
14. 1[ Indeed the argument used to show the "stay away from the cut locus" property can be 
easily adapted in this case, and all the other arguments remain true. 
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